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摘要 :目的 本 文 提出 腓 肠 肌 羽 状 角 的 一 种 基于 超声 射频 信号 的 互相 关 自 动 测量 方法 。 方 法 本 文采 集 正 常人 体 腓 肠 肌 与 仿真 
腓 肠 肌体 模 的 超声 射频 信号 ,重建 超声 图 像 , 确 定 肌 束 的 搜索 起 始点 ,采用 互相 关 算 法 自动 跟踪 肌 束 腊 和 肌肤 膜 , 绘 出 肌 束 线 和 


肌 妥 膜 线 , 计 算 两 线 之 夹 角 , 即 为 腓 肠 肌 羽 状 角 , 并 检验 自动 与 手动 两 种 测量 方法 的 可 重复 性 及 一 致 性 。 结 果 利用 互相 关 自 动 


测量 方法 可 准确 测量 出 体 膜 中 仿真 肌 束 羽 状 角 ,其 测量 值 的 均值 与 理论 值 之 间 的 差异 小 于 1"。 利 用 互相 关 自 动 测量 方法 与 手 
动 测量 方法 获得 的 正常 人 体 腓 肠 肌 羽 状 角 分 别 为 20.48°+0.47° 与 21.49°+1.79°, 两 种 方法 的 变异 系数 CV 值 均 小 于 3%, 均 方 根 
误差 均 小 于 1°*。Bland-Altman 差 值 法 分 析 结 果 表 明 两 种 方法 的 一 致 性 良好 。 结 论 本 文 所 提出 的 腓 肠 肌 羽 状 角 自动 测量 方法 


是 基于 超声 射频 信号 的 互相 关 算 法 ,其 测量 结果 与 手动 测量 结果 基本 一 致 ,在 较 严 重 的 噪声 情况 下 仍 可 跟踪 腓 肠 肌 的 肌 束 走 


向 ,自动 测量 腓 肠 肌 的 羽 状 角 。 
关键 词 : 腓 肠 肌 ; 羽 状 角 ;超声 射频 信号 ;互相 关 算法 


Automatic extraction of the pennation angle of the gastrocnemius muscles from 


ultrasound radiofrequency signals 
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Abstract: Objective We propose a cross-correlation method for automatic extraction of the pennation angle (PA) of the 
gastrocnemius (GM) muscle from ultrasound radiofrequency (RF) signals. Methods The ultrasound RF signals of the GM 
muscles in tension condition from normal subjects and the simulated ultrasound signals were collected. After the starting point 
of tracking, a fascicle was selected in the reconstructed GM ultrasound image from the RF signals, and the fascicle and deep 
aponeurosis could be automatically tracked using the cross-correlation algorithm. The lines of the fascicle and deep 
aponeurosis were then drawn and the PA was calculated. The reproducibility of the proposed method and its consistency with 
the manual measurement method were tested. Results The angles of the simulated fascicles were precisely extracted 
automatically. The difference between the experimental measurement and the theoretical values was less than 1°. The PA 
measured automatically and manually was 20.48°+0.47° and 21.49? *1.79*, respectively. The coefficient of variation (CV) of the 
two methods was less than 3% and the root-mean square error (RMSE) was less than 1°. Bland-Altman plot showed a good 
agreement between the proposed automatic method and the manual method. Conclusion The proposed cross-correlation 
automatic measurement method can detect the orientation of the fascicle and deep aponeurosis and measure the PA based on 
ultrasound RF signals with serious speckle noise. 

Key words: gastrocnemius muscle; pennation angle; ultrasound radiofrequency signals; cross-correlation algorithm 


运动 康复 训练 有 助 于 恢复 运动 功能 障碍 患者 的 中 
枢 神 经 系统 功能 ,并 可 防止 肌肉 萎缩 ", 以 减轻 运动 功能 
障碍 给 患者 及 其 家 庭 .社会 带 来 的 沉重 负担 。 运 动 康复 
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训练 时 ,不 仅 要 对 患者 施加 精确 的 运动 控制 .记录 治疗 
数据 ,更 需要 及 时 了 解 肌肉 功能 的 改善 情况 ,以 改进 训 
练 方法 和 强度 。 肌 肉 结构 特征 参数 ,如 肌肉 厚度 、 肌 纤 
维 长 度 、 羽 状 角 等 ,常用 于 评估 肌肉 功能 的 变化 1。 

超声 影像 (US) 技 术 具 有 无 创 .实时 和 便捷 的 优点 ， 
已 经 被 越 来 越 多 地 用 于 肌肉 的 研究 。 然 而 ,临床 医生 往 
往 是 通过 图 像 处 理 软件 依靠 手工 标注 的 方法 来 确定 这 
些 肌 肉 结构 参数 ,这 限制 了 超声 影像 技术 在 骨 旷 肌肉 
临床 研究 中 的 应 用 。 近 几 年 ,研究 学 者 尝试 通过 各 种 算 
法 对 肌肉 超声 图 像 进行 处 理 分 析 , 以 期 自动 获得 这 些 参 
Zi", Zhou 与 Zheng“% 提 出 迭代 霍 夫 变换 (Revoting 
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Hough Transform) 计 算 超声 图 像 中 肌纤维 的 羽 状 角 , 取 
得 了 良好 的 效果 ,并 进一步 应 用 该 算法 在 肌肉 运动 过 程 
中 动态 获取 羽 状 角 。Zhao 和 Zhang"" 采 用 雷 登 变换 
(Radon Transform) 实 现 肌肉 羽 状 角 的 自动 跟踪 和 计 
算 。 这 些 方 法 基于 超声 灰 度 图 像 ,斑点 噪声 对 跟踪 结果 
影响 较 大 ,需要 对 超声 图 像 进行 一 些 专门 的 预 处 理 , 例 
如 降 采 样 .插值 .和希 尔 们 特 变换 .幅度 -亮度 转换 KRYE 
波 等 "" ,以 获得 较 好 结果 。 而 经 预 处 理 的 超声 灰 度 图 像 
往往 会 丢失 一 部 分 超声 射频 原始 信号 所 包含 的 有 用 信 
息 ,例如 相位 信息 。 

超声 射频 信号 是 原始 信号 ,是 入 射 超声 波 经 生物 组 
织 反 射 .散射 后 形成 的 回 波 , 富 含 组 织 内 部 信息 。 与 基 
于 灰 度 网 的 跟踪 方法 相 比 ,基于 射频 信号 的 跟踪 方法 可 
利用 更 多 的 有 用 信息 , 且 受 斑点 噪声 影响 较 小 。Luo 和 
Konofagou2 提 出 一 种 基于 超声 射频 信和 号 的 快速 归 一 化 
互相 关 方 法 进行 心肌 组 织 的 位 移 估 计 , 取 得 了 较 好 的 
结果 。 

肌肉 痉挛 等 腿 部 疾病 都 与 腓 肠 肌 有 着 直接 或 者 间 
接 的 关系 。 当 肌肉 拉 伸 或 收缩 时 , 羽 状 角 会 随 之 变化 ”。 
羽 状 角 的 变化 可 反映 肌肉 的 功能 与 状态 ,为 临床 上 腓 肠 
肌 病 变 诊断 以 及 运动 康复 训练 提供 重要 的 参考 信息 。 
本 文 提出 腓 肠 肌 羽 状 角 的 一 种 基于 超声 射频 信和 号 的 自 
动 测量 方法 ,通过 选 定 肌 束 的 超声 射频 信号 ,利用 互相 
关 算法 ,分 析 射 频 信 号 的 幅 值 和 相位 信息 ,自动 跟踪 肌 


e men nM ene) 
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用 于 人 体 和 仿真 体 膜 的 超声 数据 ,所 检测 的 羽 状 角 数 值 
分 别 与 手工 测量 结果 和 理论 数值 相 比 较 ,检验 该 算法 的 
有 效 性 。 


1 材料 和 方法 
1.1 人 体 超声 数据 的 获取 

采用 开放 式 综 合 超 声 系 统 (Sonix RP, Ultrasonix, 
Canada) ,该 系统 配 有 主机 箱 ( 包 子 PIV 3.0 GHz 的 双 
CPU,80 GB 的 系统 硬盘 、 及 1GB 的 存储 器 ) 显示 器 一 
张 内 置 10 位 模 / 数 转换 卡 (采样 频率 40 MHz) ,以 及 一 
个 中 心 频率 10 MHz 的 线 阵 超声 探头 (-6dB 带宽 5~ 
14 MHz)。 该 系统 采用 开放 式 PC 平台-Microsoft® 
Visual Studio® 环 境 ,可 以 同时 保存 超声 图 像 超声 视频 
和 超声 射频 信号 。 

受 试 者 为 正常 男性 ,22~26 岁 ,坐姿 , 膝 关节 成 近 
90° , 脚 着 地 ,以 脚跟 为 定点 ,中 关节 背 届 40" ,保持 不 
动 。 腓 肠 肌 内 侧 头 为 研究 对 象 ,将 超声 探头 沿 腓 肠 肌 纵 
向 放置 。 当 显示 器 清晰 地 显示 出 腓 肠 肌 的 羽 状 肌 束 结 
构 时 ,采集 超声 射频 信号 (图 1A) ,并 重建 超声 图 像 (图 
1B)。 然 后 ,利用 本 文 所 提出 的 互相 关 自 动 测量 的 方法 
进行 羽 状 角 检 测 ( 详 见 1.3 节 ) ,检测 结果 将 与 3 位 有 经 
验 临床 超声 医师 的 手动 测量 结果 进行 比较 。 


1 腓 肠 肌 的 超声 射频 信号 、 由 射频 信号 重建 的 B- 
型 超声 图 及 仿真 腓 肠 肌 的 超声 图 
3  Fig1 Ultrasound radiofrequency (RF) signals (A), 
B-mode | image of gastrocnemius muscle 
7 4 reconstructed using RF signals (B), and a simulated 


ultrasound image of the gastrocnemius muscle (C). 


Deep aponeurosis 


1.2 仿真 体 膜 数 据 的 获取 

编写 仿真 程序 ,调用 Field I 程序 包 , 设 计 1 块 40 mmx 
10 mmx90 mm 的 腓 肠 肌体 模 , 其 中 100 000 个 散射 因 
子 随机 分 布 在 体 模 内 部 。 线 阵 超声 探头 为 128 mR, P 
心 频率 为 3 MHz, 系统 采样 频率 为 100 MHz, 然 后 利用 
Field I 计算 出 128 根 反映 射 频 信号 线 , 图 像 动态 显示 范 
围 压缩 到 40 dB。 腓 肠 肌 的 仿真 超声 结果 如 图 1C 所 示 ， 
体 模具 有 5 条 不 同 角 度 的 仿真 肌 束 和 一 条 仿真 肌 妥 


膜 。 利 用 本 文 所 提出 的 互相 关 自 动 测 量 的 方法 进行 羽 
状 角 检测 ( 详 见 1.3 节 ), 测 量 值 与 仿真 设计 时 的 理论 值 
进行 比较 ,分 析 本 文 所 提出 的 方法 的 有 效 性 。 
1.3 基于 超声 射频 信号 的 互相 关 自 动 测量 算法 

如 图 1B 所 示 , 羽 状 肌 的 肌肉 层 中 含有 浅 层 肌 妥 膜 、 
多 根 肌 束 、 和 深层 肌 腿 膜 ,根据 肌 束 的 排列 方向 和 深层 
肌 庚 膜 的 方向 可 以 画 出 肌 束 线 和 深层 肌 胖 膜 线 ,两 线 的 夹 
角 即 为 肌 束 与 深层 肌 腿 膜 之 间 的 夹 角 , 称 为 羽 状 角 "。 
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本 文 仅 以 跟踪 一 条 肌 束 为 例 介 绍 本 文 所 提出 的 羽 状 角 
测量 方法 。 

利用 Matlab7.11 自 编 “肌肉 形态 分 析 软 件 ”, 功 能 主 
要 包括 射频 信号 读 人 ,超声 图 像 重 建 .参数 选择 测量 结 
果 显 示 与 保存 等 。 首 先 , 读 人 所 采集 的 腓 肠 肌 的 超声 射 
Miei er CE CE) ,然后 ,将 超声 信号 沿 纵 向 方向 进行 降 
采样 , 沿 横向 方向 进行 插值 ,保证 重建 的 超声 图 像 与 超 
声 设备 所 显示 的 B- 型 超 图 像 大 小 一 致 (340x321) ,如 图 
1B 所 示 。 选 择 羽 状 角 参 数 , 即 可 按照 下 面 的 互相 关 算 
法 自动 测量 羽 状 角 , 并 将 结果 显示 在 结果 框 中 。 

互相 关 自 动 测量 算法 

在 超声 图 像 中 选择 一 条 清晰 的 .长 度 适中 的 肌 束 作 
为 被 测 肌 束 ,在 其 上 确定 一 跟踪 起 始点 ,以 该 点 为 中 心 
沿 该 条 信号 线 选 取 一 段 信号 ( 含 20 个 采样 点 ) 作 为 感 兴 
趣 的 肌 束 信号 x ,如 图 2A 所 示 , 然 后 沿 两 侧 相 邻 的 超 
声 射 频 信 号 线 搜索 与 x, 相关 性 最 大 的 肌 束 信号 ,搜索 
范围 设 定 为 以 该 中 心 点 两 侧 各 20 个 点 的 范围 ,如 图 2B 
所 示 。 公 式 (1) 为 互相 关 算法 公式 。 

Rrx (T) = EX X(t + 7)) (1) 

其 中 ,E 表 示 期 望 值 , x 是 当前 射频 信号 线 上 感 兴 
趣 的 肌 束 信号 , x 是 相 邻 射频 信号 线 的 肌 束 信号 ,* 表 
示 卷 积 。 搜 索 范围 为 (-20,20),t 表 示 x 相对 于 x 的 位 
移 偏 移 量 ,t e(-20, 20)。 当 R 值 最 大 时 ,表明 x, 5j x, IE 
配 , x 就 是 所 要 找 的 肌 束 信号 ,并 记录 其 位 置 [公式 
QU. 

A5-2AÀ,-TT (2) 


图 2 当前 超声 射频 信号 与 相 邻 超声 射频 信号 

Fig.2 Ultrasound RF signals of the current line and the 
neighbor line. A: Current ultrasound RF signals. x; 
presents the signal of muscle fascicle (marked by red 
box); B: Adjacent ultrasound RF signals. Two dashed 
lines indicate the searching range (-20, 20). x; presents 
the signal with the maximum cross-correlation value. 


其 中 , A, 是 相 邻 信号 线 上 肌 束 x, 的 中 心 位 , A 是 
x 的 中 心 点 的 坐标 ,t 是 x 相对 于 x1 的 位 移 偏 移 量 。 

然后 ,以 这 两 条 相 邻 超声 射频 信号 为 当前 信号 线 ， 
分 别 向 左右 方向 重复 上 述 互 相关 搜索 过 程 ,跟踪 肌 束 方 
向 。 评 判 跟踪 是 否 结束 有 两 个 标准 : 

(1) 跟 踪 到 边界 , 即 搜索 到 第 一 条 信号 线 或 者 最 后 
一 条 信号 线 ; 

(2) 跟 踪 到 肌 束 的 终点 ( 若 该 点 的 有 R 值 远 远 小 于 前 
一 个 点 的 有 R 值 , 则 该 点 为 肌 束 终点 。 本 文 设 定 肌 束 终点 
的 R 值 小 于 前 一 个 点 的 R 值 的 1/10)。 

当 跟 踪 完 成 时 ,将 各 个 信号 线 上 的 肌 束 信号 位 置 进 
行 线性 拟 合 ,得 到 跟踪 目标 一 一 腓 肠 肌 肌 束 线 。 

最 后 ,进行 深层 肌 妥 膜 线 的 跟踪 ,算法 同上 。 公 式 
(3) 和 公式 (4) 分 别 为 肌 束 线 和 深层 肌 腿 膜 线 的 线性 表 
达 公 式 , 由 斜率 值 和 便 可 计算 肌 束 线 和 深层 肌 胁 
膜 线 的 夹 角 4, 即 羽 状 角 [ 公 式 (5)]。 


y=k:x+Db (3) 
yo»-ko:x-b, (4) 
a — arctan umm (5) 


1.4 自动 手动 两 种 测量 方法 的 可 重复 性 及 一 致 性 检验 

本 文通 过 计算 变异 系数 (CV) 和 均 方 根 误 差 
(RMSE)2 个 统计 学 参数 来 评价 本 文 所 提出 的 互相 关 自 
动 测量 方法 的 可 重复 性 。 

变异 系数 是 衡量 各 测量 值 变 异 程度 的 一 个 统计 量 ， 
定义 如 公式 (6): 

CV= $x 100% (6) 


其 中 ,X 及 S 分 别 为 测量 数据 的 样本 平均 数 和 标准 
差 。 由 定义 可 以 看 出 ,CV 是 一 种 相对 分 散 度 的 衡量 指 
标 ,可 以 消除 平均 数 不 同 对 各 组 数据 变异 程度 比较 的 影 
m, 

均 方 根 误差 的 定义 如 公式 (7): 


RMSE = DN (7) 


式 中 ,n 为 测量 次 数 ,qd 为 第 i 个 测量 值 与 均值 的 差 值 。 
通过 绘制 Bland-Altman 差 值 图 ,分 析 手 动 测量 方 
法 和 互相 关 自 动 测量 方法 的 一 致 性 。 


2 结果 
2.1 仿真 实验 结果 

利用 本 文 所 提出 的 互相 关 自 动 测量 方法 检测 仿真 
超声 图 中 的 腓 肠 肌 肌 束 的 羽 状 角 ,图 3A 显示 该 算法 可 
以 很 好 地 跟踪 肌 束 及 肌 妥 膜 。 表 1 列 出 仿真 超声 图 中 5 
条 肌 束 线 与 肌 腿 膜 线 之 间 的 夹 角 ( 羽 状 角 ) 和 肌 腿 膜 角 
度 的 理论 值 和 测量 值 ,以 及 两 者 之 间 的 误差 值 。 测 量 值 
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踪 的 肌 束 线 的 长 度 往往 要 短 于 手动 绘制 结果 。 但 是 ,从 
结果 看 ,这 对 羽 状 角 的 计算 影响 不 大 ,两 种 方法 所 测量 
的 羽 状 角 平 均值 的 差 值 约 为 1" ,手动 测量 的 羽 状 角 为 
21.49°+1.79° ,自动 跟踪 的 羽 状 角 为 20.48"+0.47"。 


的 均值 与 理论 值 之 间 的 差异 在 1" 之 内 。 
2.2 临床 实验 结果 

如 图 3B-C 及 表 2 所 示 , 本 文 算法 自动 跟踪 人 体 腓 
肠 肌 的 肌 腿 膜 线 与 超声 医师 手动 绘制 结果 近似 ,而 所 跟 


Deep aponeurosis 


图 3 仿真 超声 图 和 人 体 腓 肠 肌 超 声 图 的 肌 束 线 和 深层 肌 胁 腊 线 的 检测 结果 


Fig.3 The tracking results of muscle fascicles and deep aponeurosis. A: Five muscle 


fascicles and one deep aponeurosis are automatically tracked (marked by red line) in 
the simulated ultrasound image. The five muscle fascicle pennation angles are 
marked ou-as; B: In reconstructed ultrasound image, the muscle fascicle and deep 
aponeurosis are automatically detected using the proposed method (marked by red 
lines). The pennation angle is marked a; C: In the same ultrasound image as B, the 
muscle fascicle and deep aponeurosis are drawn by hand (marked by the blue lines). 
The pennation angle is marked by a. The white curves in A and B present the traces 
detected using the proposed method. 


表 1 仿真 实验 结果 


Tab.1 Results of the simulation test 


Pennation angle 


Angle of deep aponeurosis 


Qu 02 Q5 4 Qs 

Theoretical value (°) 23.56 24.79 25.81 26.64 27.26 0 

Measurement value (°) 23.30-20.26 | 24.05x0.29 25.86+0.13 26.05+0.10  27.44x0.19 0.10+0.03 

Difference (°) -0.26 -0.74 0.05 -0.59 0.18 0.10 
A2 临床 实验 结果 
Tab.2 Results of the clinical test 

Pennation angle Angle of deep 
aponeurosis 

Manually measured value (°) 21.49+1.79 0.18+0.72 
Automatically measured value (°) 20.48+0.47 1.12+0.46 
Difference (°) -1.01 0.94 
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2.3 自动 手动 两 种 测量 方法 的 可 重复 性 及 一 致 性 检验 结果 

自动 测量 方法 与 手动 测量 方法 所 获得 的 两 组 数据 
的 变异 系数 CV(2.32% ,2.98% ) 和 均 方 根 误差 RMSE 
(0.46°,0.61°) ,分 析 两 种 方法 的 可 重复 性 。 可 看 出 两 种 
方法 的 CV 值 都 小 于 3%,RMSE 均 小 于 1°, 表 明 互 相关 
自动 测量 方法 的 可 重复 性 与 手动 测量 方法 相似 。 

图 4 显示 的 是 Bland-Altman 差 值 法 分 析 结 果 。 手 
动工 测量 和 自动 测量 两 组 数据 的 差 值 均值 与 差 值 方差 
4) 381] A -0.73? 5j 0.90? , 则 9526 — Sc PE FIR W -0.73? + 
1.96x0.90, HII (1.04° ,-2.49° ) ,其 一 致 性 界限 较 罕 , 且 只 
有 一 个 数据 点 落 到 95% 一 致 性 界限 外 ,表明 手动 测量 与 
自动 测量 这 两 种 方法 的 一 致 性 较 良 好 。 
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图 4 互相 关 自 动 测量 方法 与 手动 测量 方法 的 Bland-Altman 图 
Fig.4 The Bland-Altman plot of the pennation angle measured 
with cross-correlation automatic and manual measurements. 


3 讨论 

本 文 结 果 表 明 ,互相 关 自 动 测量 方法 可 通过 跟踪 肌 
束 的 超声 射频 信号 ,在 较 严重 的 斑点 噪声 情况 下 仍 可 准 
确 地 跟踪 腓 肠 肌 肌 束 和 肌 妥 膜 的 方向 ,自动 测量 腓 肠 肌 
的 羽 状 角 , 与 手动 测量 结果 相差 较 小 。 羽 状 角 的 大 小 易 
受 多 种 因素 影响 ,例如 受 试 者 的 个 体 差异 (年 龄 .性别 、 
体质 量 等 ) 测试 姿势 测量 位 置 . 肌 肉 的 健康 状况 、 及 超 
声 设备 的 因素 。 因 此 ,通过 限制 受 试 者 的 年 龄 体重 性 
别 ,可 减少 受 试 者 的 个 体 差 异 对 羽 状 角 测 量 的 影响 。 而 
由 于 测试 姿势 的 不 同 , 腓 肠 肌 的 状态 不 一 样 ,会 导致 测 
量 的 羽 状 角 数 值 有 所 不 同 。 本 文 的 受 试 者 取 坐 姿 , 膝 关 
5 yu yr 909 , ERA TS H JE 40 ,上 腓 肠 肌 处 于 拉 伸 状态 ， 
上 腓 肠 肌 羽 状 角 的 测量 值 为 20.48*+0.47"。 这 一 结果 与 
文献 [17] 的 测量 结果 近似 ,其 实验 对 象 采取 躁 关节 中 届 
45" 姿 势 , 健 侧 的 腓 肠 肌 羽 状 角 为 20.38"+5.42"。 最 近 有 
文献 报道 羽 状 角 的 动态 测量 ,在 路 届 过 程 中 羽 状 角 变化 
范围 在 19.3?"+0.7? 至 35.4*+6.3" 之 间 呈 。 另 有 文献 报道 
腓 肠 肌 内 侧 头 远 端的 羽 状 角 约 为 13” ,可 见 腓 肠 肌 的 
观测 位 置 对 羽 状 角 的 影响 较 大 。 


不 同 的 超声 设备 也 会 导致 与 羽 状 角 测量 结果 的 差 
异 ,主要 原因 是 成 像 算法 和 图 像 后 处 理 算法 的 差异 。 超 
声 设备 系统 在 超声 成 像 时 ,会 将 原始 超声 射频 信号 经 过 
降 采 样 或 插值 处 理 后 ,转化 为 灰 度 图 像 ,因此 不 同 的 成 
像 算法 就 会 影响 到 羽 状 角 的 测量 结果 。 成 像 处 理 后 , 往 
往 还 要 经 过 去 噪 ,滤波 等 处 理 , 以 得 到 的 清晰 的 超声 灰 
度 图 像 。 但 是 ,经 过 这 些 处 理 后 超声 射频 信和 号 所 含有 的 
生物 组 织 原始 信息 将 部 分 丢失 "”。 本 文 所 提出 的 互相 
关 自 动 测量 方法 直接 对 超声 射频 信和 号 进行 运算 ,利用 射 
频 信和 号 幅 值 和 相位 的 信息 , 受 噪 声 影响 小 。 

然而 ,这 种 基于 超声 射频 信号 的 互相 关上 自动 测量 方 
法 尚 有 不 足 之 处 ,例如 有 些 肌 束 纤维 的 连续 性 不 好 ,在 
超声 图 像 中 旦 小 段 显示 ,对 于 射频 信号 的 肌肉 纤维 , 互 
相关 自动 测量 方法 虽然 依靠 小 段 肌 束 仍 可 确定 肌 束 方 
向 ,对 羽 状 角 的 计算 影响 不 大 ,但 是 会 对 肌 束 长 度 的 测 
量 产 生 较 大 误差 ,因为 羽 状 肌 的 肌 束 长 度 定 义 为 肌 束 从 
浅 层 肌 腿 膜 到 深层 肌 用 膜 之 间 的 长 度 ""。 因 此 ,车 增加 
肌 束 长 度 测量 功能 时 ,可 基于 本 文 结果 ,利用 先 验 知识 
增 大 搜索 范围 和 改进 跟踪 终止 条 件 , 将 断 开 的 肌 束 纤维 
连接 起 来 ,然后 跟踪 计算 ,以 期 提高 肌 束 参数 测量 的 准 
确 性 。 

超声 定量 评价 对 评估 肌肉 功能 的 变化 ,以 及 对 肌肉 
运动 康复 训练 部 有 着 重要 的 意义 站 。 然 而 ,肌肉 结构 
化 参数 的 测量 目前 大 多 是 手工 标注 测量 , 面 对 大 量 的 肌 
肉 超声 图 像 , 这 样 测量 工作 费时 费力 ,自动 测量 方法 的 
研究 具有 现实 意义 。 目 前 ,所 报道 的 腓 肠 肌 羽 状 角 检 
测 方法 大 多 是 基于 超声 灰 度 图 像 ,通过 对 超声 图 像 
处 理 分 析 , 自 动 获得 肌肉 的 结构 参数 ,取得 了 良好 的 效 
果 “" ”而 本 文 直接 对 超声 射频 信号 跟踪 肌 束 方 
向 ,计算 羽 状 角 ,在 噪声 环境 中 也 得 到 了 较 好 的 结果 。 

总 之 ,本 文 提出 一 种 基于 超声 射频 信号 的 自动 追踪 
肌 束 方向 并 测量 羽 状 角 的 互相 关 算 法 ,通过 仿真 数据 和 
人 体 数据 的 检验 ,该 算法 的 可 重复 性 良好 ,与 手动 测量 
结果 基本 一 致 。 本 方法 在 较 严 重 的 噪声 情况 下 仍 可 跟 
踪 腓 肠 肌 的 肌 束 走向 ,自动 测量 腓 肠 肌 的 羽 状 角 。 在 此 
研究 的 基础 上 可 进一步 改善 以 自动 测量 肌肉 厚度 和 肌 
束 长 度 等 参数 形成 肌肉 形态 定量 分 析 软 件 , 用 于 评估 肌 
肉 功能 的 变化 ,指导 肌肉 运动 康复 的 训练 。 
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